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Abstract. Visual, NIR, and H a surface photometry of 
the oval disk galaxy M 94 (NGC 4736) was performed to 
study the distribution of mass in old stars and the distri- 
bution of the warm emission line gas. The radial profile 
of the inner stellar disk is steeper than exponential and 
displays twisted and deformed isophotes. Dust is not re- 
sponsible for these deviations. The investigation of a num- 
ber of morphological models showed that M 94 has a weak 
central stellar bar of 0.7 kpc major axis length which com- 
prises about 14% of the total light within 20". 

By means of longslit spectroscopy the kinematics of 
gas and stars in this region were investigated. The stellar 
kinematics reveal the existence of a small spheroidal bulge 
with v/a 1=3 0.8. The stellar velocity field in this region 
is approximately axisymmctric, showing that the effects 
of the bar on the kinematics of the hot component are 
relatively small. The warm gas in the bar region shows 
global and local deviations from the stellar kinematics, 
but outside the central bar fits into the general HI velocity 
field. 

Model calculations of closed orbits for gas in the po- 
tential of the bar perturbation and axisymmetric disk and 
bulge predict large non-circular motions for cold gas in 
equilibrium flow. These do not fit the central H a kine- 
matics; rather it appears that the H a gas must have large 
random motions and is not in a steady state, and that hy- 
drodynamical effects due to a recent starburst play a role 
in the central region. 



1. Introduction 

M 94 = NGC 4736 is an Sab spiral galaxy at a distance of 
« 6.6 Mpc. It is well studied in optical morphology, optical 
emission gas kinematics (van der Kruit 1974, 1976, Buta 
1984, 1988) and HI distribution and kinematics (Bosma 
et al 1977, Mulder & van Driel 1993). There exist several 
good overviews about these results and other characteris- 
tics of M 94 in the literature (e.g. Buta 1984, Mulder & 
van Driel 1993). Since this paper deals mainly with new 
optical and NIR results (morphology, stellar kinematics, 
gas kinematics) we present here only a brief summary of 
the most important optical features of this galaxy. Basic 
parameters of M 94 and the most numerical results from 
this study are presented in Table 1. 

Bosma et al (1977) distinguished the following five an- 
nular zones with progressivly lower surface brightness in 
M94: 



A bright central region (bulge), R < 15". 
A zone of inner spiral structure R ~ 15 to 50". This 
zone is bounded by a bright inner ring. 
A zone of outer spiral structure, R w 50 to » 200", 
this region probably forms an oval disk. 
A zone of low surface brightness (gap). 
A faint outer ring of R rs 330". 



Key words: 



Send offprint requests to: C. Mollenhoff 

* Visiting astronomer of the German-Spanish Astronomi- 
cal Center, Calar Alto, operated by the Max-Planck-Institut 
fur Astronomie, Heidelberg jointly with the Spanish National 
Commission for Astronomy. 



The nuclear region is extremely bright; according to 
Keel and Weedman (1978) M 94 is one of the galaxies with 
the highest optical surface brightness. The inner ring of 
radius rj 40"to 50", shows HII regions and young blue 
objects, characteristic for strong star formation. The in- 
terpretation of this ring is the subject of numerous papers. 
In a detailed study of the kinematics van der Kruit (1974, 
1976) attempted to fit a variety of kinematical models to 
the ring. He concluded that the velocity field is best inter- 
preted in terms of a uniform radial expansion from the nu- 
cleus. However, using Fabry-Perot data Buta (1984, 1988) 
concluded that the ring is better interpreted as the result 
of a secular evolution in this galaxy due to a Lindblad 
resonance. 
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The bright inner ring was also observed in CO by Gerin 
et al (1991). They followed the formation of the outer and 
the inner ring in a model where the main disk of M 94 is a 
massive oval with b/a — 0.8. That the disk of M94 is not 
axially symmetric had already been seen by Bosma et al 
(1977). 

The zone inside the inner ring displays a strong 
isophote twist. This could be the consequence of a tri- 
axial bulge (Beckman et al 1991) or a central secondary 
bar (Kormendy 1993). Such a 'bar within a bar' would ro- 
tate much faster than the large outer oval disk (Shlosman 
et al 1989, Friedli and Martinet 1993). A third interpreta- 
tion was given by Shaw et al (1993), who proposed that 
this structure is the result of a combined stellar- and gas- 
dynamical mechanism. In the potential of the rotating oval 
disk a gas ring formas which perturbs the central stellar 
distribution. This pertubation is misaligned with respect 
to the outer oval disk but has the same angular speed; 
thus it is not an independent bar. 

The aim of this paper is to study the distribution 
and kinematics of stars and gas in the central region of 
NGC4736, and to compare these with model calculations. 
In Sect. U we describe the observations: (1) Surface pho- 
tometry in the visual region (V,I,H a ) and in the NIR 
(J,K), in order to study the distribution the old stellar 
component as well as that of the dust and warm, line- 
emitting gas. (2) Longslit spectroscopy for the investiga- 
tion of the kinematics of stars and gas in this region. The 
data reduction and the results are described In Sect. |[ 
In Sect. H we present morphological models for M 94. In 
Sect. ^ orbit families in a corresponding model potential 
are calculated and are used to model the kinematics of 
stars and gas. Their predictictions will be compared with 
the observed kinematics. In Section ^| we discuss our re- 
sults, and the conclusions of this paper are given in Sec- 
tion . 

2. Observations 

2.1. CCD surface photometry 

The CCD imaging observations were performed in Febru- 
ary 1990, in January 1992, and in July 1993 at the 
German-Spanish Astronomical Center on Calar Alto, 
Spain. The Calar Alto 1.23m telescope was used in three 
different setup modi for CCD imaging: (1) The Lan- 
dessternwarte focal reducer at the Cassegrain focus with a 
GEC P8603/Achip (576x386 pixels, 1 pix of 22^ = 1.56") 
and V, I filters; (2) the Calar Alto Cassegrain focus CCD 
camera with a coated GEC chip (576 x 386 pixels, 1 pix of 
22^ = 0.46") and / filter; (3) the Cassegrain focus CCD 
camera with a Tektronix chip (512 x 512 pixels, 1 pix 
of 27/z = 0.56") and H a filter. The exposure times were 
1 minute (/ filter) and 3 minutes (V) with the focal re- 
ducer, 15 minutes (/) and 45 minutes (H a filter) with 
the Cassegrain camera. Photometric standard stars were 



exposed several times during the night for the flux cal- 
ibration. Flat field exposures were taken during evening 
and morning dusk; the telescope was slightly displaced 
between each exposure, to allow the elimination of stellar 
images in the flat field frames by an appropriate filtering 
procedure. 

2.2. NIR surface photometry 

The NIR imaging observations were performed in January 
1994 with the 2.2m telescope of the German-Spanish As- 
tronomical Center on Calar Alto, Spain. J, K images were 
obtained with the MPIA MAGIC camera (Herbst et al 
1993). MAGIC has a 256x256 NICMOS3 HgCdTe detec- 
tor array and was used in an imaging mode of 0.66"/pix, 
corresponding to a field of 169"xl69". 3 exposures of 20 
sec were added in the camera electronics, using a 'sam- 
ple up the ramp' readout mechanism. Then the telescope 
was moved by a few arcsec for the next exposure, and 
another few arcsec for a third exposure, in order to get 
rid of the bad pixels. In total six such triples, i.e. 18 im- 
ages of 1 min exposure time were obtained in each filter. 
Between the galaxy exposures the telescope was moved 
approx. half a degree offset for 3 blank sky images of the 
same expore times, also shifted slightly from each other 
to remove background stars. Photometric standard stars 
were exposed for the flux calibration, the sky exposures 
serve as flat fields. 

2.3. Kinematics of the stars 

The spectroscopic observations were done during Febru- 
ary 1990 with the 2.2m telescope of the German-Spanish 
Astronomical Center on Calar Alto, Spain. Longslit spec- 
trograms were obtained using the standard Boiler & 
Chivens spectrograph in 4 slit orientations, P. A. = 
0°, 45°, 90°, 135° respectively. The detector was an RCA 
CCD with 1024 x 640 pixels of 15/i size. The slit width 
was 250/i or 2.9". 

For the stellar kinematics we used a grid with 1200 
lines/mm in the wavelength range 4500 - 5400 A reso- 
lution 0.9 A/pix). 3 pixels (or 2.7") along the slit were 
binned during readout. The exposure times were 2 x Ih 
for each slit position. 

Several stellar templates of spectral type between G8 
III and K5 III were exposed for the subsequent Fourier- 
Cross-Correlation. The internal He-Ar lamp was used for 
wavelength calibration. A large number of domeflats with 
different exposure times were exposed for the correction of 
the nonlinear columns of the RCA chip at low light levels. 
Skyflats were obtained to correct for the profile along the 
slit. 

2.4- Kinematics of the gas 

The spectrograms for the kinematics of the emission gas of 
the HII regions were obtained during the same observing 
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run with the same slit width and orientations as above. 
The grid for the gas kinematics had 1200 lines/mm in the 
wavelength range 6200 - 7100 A, and the exposure times 
here were 30 min for each spectrogram. The full spatial 
resolution of 0.9 "/pix along the slit was used. 

During a second spectroscopic run in January 1994 
gas spectrograms in further slit orientations P. A. = 
28°, 73°, 118°, 163° respectively, were obtained with the 
same spectrograph and grating at the same telescope. 
However, this time a Tektronix chip 1024 x 1024 of 24/i 
pixelsize, slit width 370/^ or 4.3", and spatial resolution 
of 1.47"/pix was used. The above given slit orientations 
correspond to the position angle of the central bar in M 94 



(Sect. 3.1) 



3. Data reduction and results 

3.1. CCD surface photometry 

The CCD images were reduced using the standard 
procedure with bias subtraction and flatfield division. 
An intensity-dependant median filter was applied which 
smoothed the images only in the outer zones of low S/N 
ratio but not in the central region of the galaxy. 

Figure |l| shows the central 2x2 arcmin of the / image 
of M 94. The slightly irregular outer structure is the well 
known oval ring of HII regions, not very pronounced in 
this deep red image. While the major axes of the outer 
isophotes at r w 40" have a position angle (P.A.) w 120°, 
they are nearly circular at r ~ 20". Further inside the 
position angle changes to w 28°, and the isophotes dis- 
play a distinctly cusped structure. Beckman et al (1991) 
interpreted the isophote twist as a signature for a triaxial 
bulge. However, the occurrence of cusped isophotes makes 
it more plausible that a bar-like structure of ~ 30" = 1.0 
kpc length is superposed onto the inner disk of M 94. This 
bar was already mentioned in Kormendy (1993). 

A third interpretation was given by Shaw et al (1993), 
who proposed that this structure is the result of a com- 
bined stellar- and gas-dynamical mechanism: The oval 
disk leads to an accumulation of gas in a non-axisymmetric 
circumnuclear ring which becomes phase-shifted with re- 
spect to the original stellar oval disk. This ring has suffi- 
cient mass to induce a triaxial perturbation in the central 
stellar component. This perturbation is misaligned with 
respect to the outer oval disk but has the same angular 
speed; thus it is not an independent bar. 

How does the isophote twist model of Shaw et al 
(1993) compare with our observations? Shaw et al give 
only rather general statements about the shape of the 
isophotes to be expected; it is difficult to make a com- 
parison with the observed morphological details. They in- 
terprete the region of isophote shift as a perturbation due 
to the (phase shifted) gas ring; thus it should extend out 
to the radius of the ring. However, the most important 
point is that the bar ends at a finite radius inside the ring 



(cusped isophotes). Therefore we consider the interpreta- 
tion of Shaw et al as less plausible. 

In Sect. H we will present models of the central pho- 
tometry and discuss the nature of the barred component 
in more detail. We will give arguments that it is indeed 
a bar. In the following we will therefore call it a 'bar' for 
(semantic) simplicity. 

The I image was calibrated using the standard stars 
(M 67) and comparing with the aperture photometry from 
de Vaucouleurs & Longo (1988). 



Fig. 1. Central region of M94 (NGC4736), I filter, 78 x 78 
arcsec, north is up, east to the left. A logarithmic saw-tooth 
lookup table was used to produce this image, which shows iso- 
contours as dark lines and the shallow gradients in between 
as different gray shadings. The cusped deformations of the 
isophotes due to the central bar at P.A. = 28° are easily visible. 



3.2. NIR surface photometry 

The data reduction of the J, K images proceeded as fol- 
lows: To each triple of sky exposures a vertical median fil- 
ter was applied, thus removing the stars. These sky frames 
were obtained for each galaxy and each filter, they were 
subtracted from the 18 galaxy images respectively. Since 
the dark current of the MAGIC camera is very low (Herbst 
et al 1993) it was not corrected for. The internal accuracy 
between all sky exposures of one night in each filter is bet- 
ter than 0.2%, with larger deviation from night to night. 
Therefore the flat field images were constructed for ev- 
ery night separately, using all sky images in each filter of 
that night. After flat field division of the 18 galaxy frames 
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the relative offsets of the galaxy centers and of stars (if 
present) were measured. Then the images were added to- 
gether pixelwise, applying their individual offsets. An in- 
ternal comparison allowed us to reject the bad pixel values. 
The J, K images look very similar to the I image (Fig. |]J), 
showing the same behaviour of the isophotes and the bar 
feature. This is another support for the statement that the 
bar is not a feature due to dust absorbtion. 



3.3. Isophote fits 

After masking of disturbing field stars the isophotes of 
the /, J, and K images were fitted by ellipses (Bender & 
Mollenhoff, 1987). The results for the K image and the 
/ focal reducer image are shown in Fig. || in the left and 
right panels respectively. The upper plots show the surface 
brightness in the inner part (left panel, K, r < 90"), and 
out to the largest radius we could obtain (right panel, /, 
r < 450"). As already stated by Boroson (1981) and Kor- 
mendy (1982), M94 shows a steep non-exponential profile 
in the inner disk region (r < 200"). An exponential pro- 
file can be fitted only in the very outer disk. However, 
this region of steeper profile is still a disk in the kinemat- 
ical sense: as will be seen later from the stellar velocity 
dispersion data, the real bulge of M94 extends only to 
ps 20" radius. The steep disk profile is not totally smooth, 
it shows two 'shelves' at r ps 45" and r ps 120" which, 
according to Kormendy (1982) are characteristics of oval 
disks. These 'shelves' can be detected in all our visual and 
NIR surface brightness profiles. 

The next plots in Fig. || show ellipticities e = 1 — b/a 
and position angles of the isophote fits to the K and / 
images. The curves are very similar in the different colors 
which means that dust does not play a role in M94 (cf. 
Sect. |3.6| ). Due to the bar structure the ellipticity profile 
has a first maximum of e = 1 — b/a = 0.22 at ps 10" radius 
(Fig. ^, left panel), the position angle is P. A. = 28° there. 
In the transition zone between central bar and outer disk 
at r ?s 20" radius the isophotes are nearly round (e ~ 0). 
Then e increases again, accompanied by a strong variation 
of P. A. between 20" and 30". With increasing radius, the 
ellipticity increases generally from 0.2 to 0.25 at r ps 120", 
and to ~ 0.4 at large radii. The general slope of 1 — b/a 
is only interupted by two local maxima in the ellipticity 
at 50" and 125", followed by a slight local decrease of e. 
These maxima correspond to the two shelves of the surface 
brightness at 45" and 120". The main point here is, that 
e ps 0.18 between r = 50 and 100" and e ps 0.23 between 
r = 120 and 220". As Kormendy (1982) mentioned, these 
nested elliptical regions of different e indicate, that at most 
one of these zones can be axial symmetric. Since warps can 
be excluded at this high surface brightness (Fig. §), part 
of the disk must be oval. If we assume the zone between 
r = 50 and 100" as axial symmetric, we get an incliation 
angle of i = 35°. Then the outer zone between r = 120 



and 220" must be intrinsically elliptic with an axial ratio 
of b/a w 0.9. 

The bottom plots of Fig. || show the a4-Fourier co- 
efficent which quantifies symmetric deviations of the real 
isophotes from ellipses (Bender & Mollenhoff 1987). A pos- 
itive 04 value marks cusped isophote deviations; they are 
especially strong at r = 10" (left panel, K). 04 is large and 
positive in the limited radius between r — 5" and 15"; this 
suggests that the bar ends at r ps 15", well inside the star 
formation ring at r ps 45". Then 04 shows an oscillating 
structure with increasing radius; at 45" and 110" we ob- 
serve a change from slightly negative 04 values (i.e. boxy 
deviations) to positive values. This finding is another con- 
fimation of the structural changes in the disk of M 94 at 
r ps 45" and r ps 120" (right panel, I). 

3.4- Fit of the surface brightness profile 

To study the light distribution in the inner part of M 94 we 
used the I Cassegrain focus image, since with 0.46"/pix 
it has a better pixel-resolution than the J, K images with 
0.66"/pix. However, the / image is still smoothed by atmo- 
spherical seeing, which was determined from stellar images 
to a value of 1.76", 

In first test fits of the surface brightness a de Vau- 
couleurs profile was found too steep while a Hernquist 
(1990) profile already gave a fairly good fit. For a better 
modelling we proceeded in a more general way, making use 
of the 7-profiles described by Dehnen (1993). They are 
characterized by a central luminosity density slope ~ r~ 7 
and a length scale a. They include profiles between a steep 
Jaffe type (7 = 2), de Vaucouleurs type (7 = 1.5), Hern- 
quist type (7 = 1), to a flat core type profile (7 = 0.5). 

For a number of different values of slopes 7 and length 
scales a we constructed a series of two-dimensional im- 
ages from a model disk of the corresponding profile, seen 
under the inclination of M 94 (i = 35° determined from 



the surface photometry of the inner disk, see Sect. 3.3). 
For this modelling we chosed a high pixel-resolution of 
1 pix = 0.25". These model images were then convolved 
with a Gaussian seeing profile of FWHM = 1.76". Then 
an ellipse fit was applied to every model image. In this 
way we obtained a two-parametric manifold of curves 
(7 = 0.7,0.8, ...1.7 and a = 12,15,17,20,25,30") which 
we compared with the observed ellipse fit curve. Fig. || 
shows the profiles of a 7-series for a — 15" against log r 
compared to the observation. The best fit for the inner 
region r < 30" — 40" was the curve for 7 = 0.9 ± 0.05 
and for a = 15 ± 1.5". This is very similar to Hernquist 's 
(1990) luminosity profile, which corresponds to 7 = 1. For 
the outer part of the model disk we added a simple expo- 
nential profile with r e = 60" and neglected the oval shape 
of the outer disk. This is justified since our model disk 
is intended mainly for the inner region of M 94, near and 
inside the inner ring. 
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Fig. 2. Results of the isophote fits to the K image (left panel) and the focal reducer / image (right panel). The plots show 
surface brightness, ellipticity, position angle, and the a4-Fourier coefficient over radius. The central bar causes a first maximum 
of the ellipticity at w 10" radius with a position angle P. A. — 28° (left panel, K). For r > 45" (outer disk) the ellipticity 
increases slowly from 0.2 to 0.4, the P. A. decreases from 120° to 90° (right panel, /). The structural changes in the disk at 
~ 45" and ~ 120" can be seen in all plots. 
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Fig. 3. Seeing-convolved 7-profiles with a = 15" for different 
7-values compared to the observed /-profile. The curve for 
7 = 0.9 gives the best fit for r < 40". 



3.5. Distibution of the optically emitting gas 

What is the distribution of the warm ionized gas in M 94? 
For that purpose a CCD image with an H a filter was ex- 
posed. Since this image still contains a strong continuum 
contribution, we used the I frame as our most appropri- 
ate continuum image. Since both images have different 
pixel sizes, we rebinned the H a image (27// pixels) to the 
pixel size of the I image (22/x). Then both images were 
centered to ~ 0.1 pix and were subtracted with an appro- 
priate intensity scaling which made the pure continuum 
regions vanish. Fig. |] shows the result, it is very similar 
to Fig. 4d in Pogge (1989). The main H a contribution 
comes from the well known ring of star forming HII re- 
gions. There exists also a central component of irregular, 
elongated shape; this extends mainly parallel to the stel- 
lar bar (compare also the H a map in Duric & Dittmar, 
1988). A weak and also irregular spiral structure can be 
detected between the ends of the bar and the strong ring 
at r = 45". The defects in the very center are due to the 
different seeing during H a and continuum exposures. 

3.6. Color index images 

V — I, I—K, and J — K, color index images were produced 
to study the morphologhy of dust absorption in M 94. The 
two respective images were rebinned, if necessary, to the 
same pixel scale. Then they were adjusted in position to 
an accuracy of < 0.1 pixels and divided in order to obtain 
a color index image. Dust reddening is present in the dark 
shaded regions in the color images. Fig. || shows, as exam- 
ple, the I—K image. There is some dust present in circular 
arcs around the center, especially in the NE and W. These 
features seem to be correlated with the faint spiral arms 



Fig. 4. H a image of M 94, 2x2 arcmin, north is up. The contin- 
uum was subtracted using an appropriate red image. Emission 
of warm gas is present in the ring and in a central region which 
coincides roughly with the bar. 

around the central bulge. They can also be detected in 
the B — V, V — R images in Beckman et al (1991, Fig. 
7a, b). Further out two spiral-arm- like dust features follow 
the ring of HII regions. All these dust structures can al- 
ready be seen in the B broad band image in Beckman et 
al (1991, Fig. 8). The most important point here is that 
the inner dust arcs are situated near the long axis of the 
bar; they would not be able to produce the bar feature 
by an adequate extinction of a non-barred central stellar 
distribution in M 94. 

3. 7. Stellar kinematics 

The CCD images of the stellar absorption spectra (4500- 
5400 A) were reduced using the standard procedure with 
bias subtraction and flatfield division. A special mask 
function was constructed from a series of many differently 
exposed dome flats in order to correct for the intensity- 
dependent additive column offsets of the RCA chip at low 
light levels. The light distribution along the slit was cor- 
rected to a flat profile by means of an apropiate func- 
tion obtained from skyflat exposures. The images were fil- 
tered to remove cosmics and then rebinned row by row in 
log A (equidistant in Av). The spectra were centered and 
placed exactly parallel to the rows. During observation 
we had obtained two spectrograms for each slit position. 
Now the reduced spectrograms of each pair were coad- 
ded with an algorithm which removed further defects and 
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Fig. 5. Color index image I — K of the central parts of 
NGC4736 (M94), 2x2 arcmin, north is up. There exist 
some dust-reddened arc regions near the center (dark shaded), 
however they are not responsible for the deformations of the 
isophotes in the I image. Long dust lanes of spiral arm shape 
are crossing the ring of HII regions at r ftf 45". 



remaining cosmics. The continuum was subtracted by ap- 
plying a polynomial fit to every row. The Hp- and [OIII]— 
lines were masked to avoid disturbances by the emission 
lines. Then the spectrograms were prepared for the cross- 
correlation by applying a cos-bell funtion to their ends. A 
similar technique was used for the preparation of the tem- 
plate spectra of several G8III, KOIII, K3III, and K5III 
stars. 

The stellar rotation velocity and the velocity disper- 
sion of the galaxy were determined by a Fourier-cross- 
correlation technique with a standard template star; it is 
described in detail in Bender (1990). This technique allows 
one also to obtain the line-of-sight velocity distribution of 
the stars. 

Fig. ^ shows the resulting stellar rotation (lower panel) 
and stellar velocity dispersion (upper panel) along the four 
slit orientations P.A. = 0°, 45°, 90°, 135°. The strongest 
rotation amplitudes are seen along P.A. = 90° and 135°. 
Since the latter amplitudes are nearly equal, the kinematic 
line of nodes must be at ~ 113° and the rotation axis at 
P.A. « 23°. The stellar rotation velocity increases nearly 
linearly to r = 12"and then remains on that level until 
r ~ 60". For larger radii our signal is too weak for a reli- 
able result from the cross-correlation analysis. The curves 
become noisy for r > 40", this is due to template mis- 
matching in the ringlike zone of star formation. 



Under the assumption of circular symmetry the stellar 
velocity field can be easily deprojected. If this assump- 
tion is justified, and all 4 rotation curves are deprojected 
with the same correct inclination with respect to the cor- 
rect line of nodes, than the deprojected rotation veloc- 
ity curves should fall together. If the adopted inclination 
and line of nodes are chooesen incorrectly, the projected 
curves are different. In this way we determined by itera- 
tion i = 30° ± 5° for the inclination, and a position angle 
of P.A. = 113° ± 1° for the line of nodes. The deprojected 
rotation curves coincide quite well, especially outside the 
bar. From the superposed projected curves it is possible 
to calculate a mean stellar rotation curve with a better 
signal-to- noise-ratio (see next Section). This 'kinematic' 
inclination angle is very similar to that obtained by sur- 
face photometry of the inner region of the disk [i = 35° 
between r = 50 and 100", see Sect. 3.3). 

In the stellar velocity dispersion profile (upper panels 
of Fig. |6|) we see a rather sudden increase inside r = 20" 
to the central value of a — 120 kms~ . This is most plau- 
sibly interpreted as the signature of a small central bulge 
with r « 20" in M 94. Outside the bulge we obtain a = 60 
kms -1 , a value which is determined by our spectral reso- 
lution and is probably higher than the velocity dispersion 
of the disk. The extent of the bulge is fairly independent 
of the slit orientation, therefore its shape is nearly axially 
symmetric. 

The broadening function of the stellar absorbtion lines 
in a galaxy spectrum is produced by the local velocity dis- 
tribution integrated along the line of sight. The broaden- 
ing function is gaussian only for simple dynamic systems; 
in general it is a more complex function, which we call the 
line-of-sight stellar velocity profile (LOS-SVP). we use the 
method of Bender (1990) to determine this function by a 
deconvolution of the correlation peaks. Fig. ^ shows the 
line-of-sight stellar velocity profiles of the spectra along 
P.A. — 90° at different radial distances from the center 
(0", ±5", ±10", ±20"). While the central profile is broad 
and gaussian (bulge profile), the line profiles at ±5" are al- 
ready asymmetric. At ±10" the profiles are more complex: 
we see a narrow peak from the rapidly rotating disk-stars, 
superposed onto a broader, more slowly rotating compo- 
nent from the bulge stars. At this radius the contributions 
of disk and bulge are approximately of the same magni- 
tude. At ±20" only the disk profiles are visible on both 
sides. 

These line-of-sight velocity profiles confirm the exis- 
tence of a bulge of ~ 20" radius in M 94. Kormendy (1993) 
used the old kinematical data of Pellet & Simien (1982) 
to argue from v ro t/cr w 0.75 and e w 0.1 that M94 has no 
bulge but is disklike even in the center. This argument is 
no longer valid since we now know that it compares ve- 
locities from two different stellar components. From the 
velocity decomposition in Fig. [?] we see that at r = ±10" 
the bulge component has a much smaller rotation veloc- 
ity than the disk component; thus v ro t/u for the bulge 
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Fig. 6. M94: Stellar velocity dispersion (top) and stellar rotational velocity (bottom) along P. A. = 0°, 45°, 90°, 135°. The 
kinematical influence of the small bulge extends only to r ~ 20" = 600pc. 



stars becomes much smaller, and is not in conflict with 
the classical v rot /<J — e - diagram. Furthermore it should 
be emphasized that spiral structure is not detectable in- 
side the bulge radius (20") of M94 (e.g. Fig. |). 

3.8. Emission gas kinematics 

The reduction of the emission line spectra was performed 
in a very similar way as described above for the stellar 
spectra. In order to extract the emission lines from the 
stellar continuum, the latter was subtracted via a poly- 
nomial fit (with masked emission lines) in every row. The 
velocity field of the gas is then easily measurable in the 
H a , [Nil], and [SI I] emission lines. The radial velocity of 
the gas was determined by cross-correlation with the emis- 
sion line spectrum from the galaxy center. Since there the 
H a emission is compensated by the absorbtion of the stel- 
lar continua, we also used as templates spectrograms from 



the star forming rings (and correcting for their rotational 
offset). 

Fig. H shows the resulting gas rotation curves for all 
8 observed slit orientations. To compare the gas veloci- 
ties with the stellar rotation field we use the mean stellar 
rotation curve from the last Section and project it along 
the corresponding slit position angles. This results in the 
solid lines in each panel of Fig. ||. The rotation curves of 
gas and stars show a similar general trend, i.e. practically 
disk dominated rotation for both components, however, 
with some marked differences. They are best studied by 
subtraction of the mean stellar rotation from the gas ve- 
locity curves for each slit orientation. 

(1) Inspection of the velocity residuals shows that 
the gas velocities are systematically higher than the stel- 
lar mean velocities. These systematic differences are not 
an artefact of our mean stellar rotation curve, they can 
be seen directly by comparing measurements of gas and 
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Fig. 8. M 94: Rotational velocity of the gas (circles) along 8 different slit positions, marked in each panel. The bulk rotation of 
the gas is similar to the stars (solid line) . However, the gas rotates faster and shows additional humps and dips in the curves. 
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Fig. 7. Line-of-sight stellar velocity profiles (LOS-SVP) along 
P. A. = 90° at different radial distances from the center (0", 
±5", ±10", ±20"). While the central profile is broad and gaus- 
sian (bulge profile), the influence of the bulge decreases with 
distance from the center. The combined profile at ±10" shows 
the broad, non-shifted bulge profile and the rotationally shifted 
narrow disk profile. At ±20" only the narrow disk profile can 
be seen. 



stellar rotation for the common position angles P. A. = 
0°,45°,90°, and 135°. Several interpretations are possi- 
ble. 

(i) The gas rotates intrinsically faster than the stars. 

(ii) A second possibility could be a different rotation 
axis for stars and gas. As we have seen in the last section, 
the best-fitting inclination of the stellar rotation plane is 
30° ±5°. A simple increase of the inclination of the gas or- 
bits to w 45° and a slight change of the line of nodes could 
explain the general trend (assuming circular gas orbits), 
but this would mean that the gas rotates in an inclined 
orbit with respect to the stellar rotation which is a rather 
implausible suggestion. 

(iii) A third possibility is that the gas moves on elliptic 
orbits along an axis roughly parallel to the small central 
bar. Deviation of the gas orbits from circular motion were 
already discussed by van der Kruit (1976), when he mea- 
sured the gas velocities in the H a -i'mg. Strong deviations 
from the pure circular motion of the HI gas were also ob- 
served by Mulder & van Driel (1993) in the form of varia- 
tions in inclination and the line of nodes of the HI velocity 
field. Our observed maximal amplitude of 120 km/sec for 
the gas velocities at P. A. = 118° between r = 20" and 40" 



matches well to the corresponding HI velocities in Mulder 
& van Driel (they use i = 40° in their Fig. 12b). 

A further plausible explanation of the systematic dif- 
ferences between stellar and gas rotation is given in Sect. 
5. 

(2) In the gas rotation curves there appear 'overshoot- 
ing' and 'counter-rotating' humps and dips of the gas ve- 
locity near the center of the galaxy. Fig. [)] shows the loca- 
tion of these local deviations by filled (humps) and open 
(dips) circles. Since many of them occur at similar radii 
(10" to 20") for different position angles, we conclude that 
the reasons for these local deviations are of nearly circu- 
lar shapes, probably spiral arms. The eastern arc of filled 
cirles (redshifted deviations with respect to the general 
velocity field) can roughly be identified with the corre- 
sponding dust arcs displayed in the I — K image (Fig. ||). 
However, there also seems to be some connection with the 
bar: The deviations occur mainly along the long sides of 
the bar and they are mostly positve on the eastern side 
and negative on the western side. In contrast to the gen- 
eral trend these deviations can not be eliminated by a 
varition of the inclination angle and/or the line of nodes. 

A comparison with the dynamical model will be given 
in Sect. 5J:, the discussion of the results in Sect. ||. 



4. Morphological model: Central bar or triaxial 
bulge? 



In Sect. 3.1 we described the isophote twist, the variations 
of the cllipticity, and the occurrence of cusped isophotes in 
M 94 inside 40" (Fig. |l] and Fig. |). It could be shown from 
the I — K color index image that these are not an effect of 
dust absorption. Are the isophote twists etc. the signature 
of a bar which is a disk phenomenon, are they caused by 
a triaxial bulge, i.e. a structure extended in z-direction? 

For the construction of a morphological model of the 
inner region of M 94 we have to consider the following 
observational facts: 

a) isophote twist of nearly 90° between 15"and 25" 

b) round isophotes inside 10"which can not be ex- 
plained just by seeing 

c) cusped isophote shapes atrw 15" 

d) increased stellar velocity dispersion inside 20", in- 
dependent of the position angle. 

We present two different models for the interpreta- 
tion of these observational findings: (1) M 94 consists of a 
disk and an axisymmetric (spheroidal) bulge, the isophote 
twist and the cusped shapes of the isophotes are due to 
a weak additional bar component. (2) M94 has a triaxial 
bulge in the center which is responsible for the isophote 
twist (no direct explanation for the cusped isophote devi- 
ations). 
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tometry of the inner disk region (Sect. 3.3): inclination 
angle i — 35°, position angle P. A. = 113°. 

This model was convolved with a Gaussian of FWHM 
= 1.76" according to the seeing during observation (mea- 
sured from the point spread function of stars in the I- 
image) . For the outer disk we adopted an exponential pro- 
file with r e = 60", with inclination and position angle as 
above. The central intensity was adjusted to the I image. 

When this model was subtracted from the observed I 
image, a central barred component remained, with a fairly 
flat profile. Figure |l(] shows the residual central bar. The 
size of the bar is about 30 x 14" in projection. At the 
distance of NGC4736 (D = 6.6 Mpc, 1" = 32pc) and for 
an inclination of i — 30° , the size of the bar from end to 
end is 42 x 14" or 1.3 x 0.45 kpc. 



Fig. 9. M 94: Locations of the humps and dips in the gas rota- 
tional velocities relative to the mean stellar rotational velocities 
along 8 different slit positions. Open circles refer to negative, 
filled circles to positive residual velocities, their diameter gives 
the absolute values. The deviations concentrate along circular 
arcs, probably connected to the inner spiral structure. 



4--1- Axially symmetric bulge, disk, and bar 

This model with an axisymmetric bulge is supported by 
the fact that the stellar velocity dispersion profile is iden- 
tical for 3 of 4 slit orientations. The a - profile along 
P. A. — 135° is somewhat broader; however, we do not 
consider this as significant. 

Since we are mainly interested in the structure of the 
inner kpc, we assume that the entire disk of M 94 is ax- 
isymmetric and neglect its outer oval shape. The mod- 
elling proceeded as follows: According to the isophote el- 
lipse fits a seeing-convolved 7 = 0.9 (Dehnen & Gerhard 
1994) model for the inner disk (plus bulge), and an outer 
exponential disk component were constructed. This two- 
component model together with an additional constant 
value for the sky background fitted the whole disk of M 94 
very well. The model parameters for the inner disk were: 
7 = 0.9, length scale a — 15", and from the surface pho- 



Fig. 10. Residual central bar component in M94 after the 
subtraction of a two-component model for bulge (7 = 0.9) and 
exponential disk. The size of the bar is 30 x 14" in projection. 
The scale is the same as the I image in Fig. |l| (78 x 78", north 
is up. 

In a second step the central bar itself was modelled. 
As Fig. [n] shows, the luminosity profile along the bar is 
rather shallow in the center and steep at the boundaries. 
Therefore an elliptical r 1 / 4 profile is not adequate. A flat- 
ter Ferrers profile (Ferrers 1877) was chosen. It has the 
following elliptical surface brightness distribution: 
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where 

The parameter n determines the radial density profile, So 
is the central surface brightness. The long and short axes 
of the bar are denoted by a and 6, respectively. For our 
model we chose n = 2, a = 15" and b = 7" (in projection). 

Fig. [ll] shows the coadded model disk and Ferrers bar. 
Subtraction of this composite model from the / image left 
only very small residuals. 



Fig. 11. Two-component model for the disk and bulge together 
with a Ferrers model for the central bar component in M94. 
Compare this composite model with the I image in Fig. jjj; the 
same scale and lookup table were used. 



4-2. Triaxial bulge and disk 

One of the main characteristics of the above model is the 
fact that the overwhelming portion of the mass in the 
central region is contained in the axisymmetric compo- 
nent, which results in a dominant monopole component 
of the potential in the disk plane. Compared to this the 
quadrupole component due to the bar is only a minor per- 
turbation of the order of a few percent. If the light in the 
inner 10" is dominated by a single triaxial component, this 
can therefore only have a small flattening in the disk plane 
of M 94. 

We have the following observational boundary condi- 
tions for the construction of a triaxial bulge: 

a) The long axis of the bulge is along P. A. = 28°. 



b) From the line-of-sight-velocity decomposition 
(Fig. 0) we conclude that the amount of light from the 
bulge at r w 10" is approximately equal to that of the 
disk there. 

c) From the a - profile (Fig. ||) the bulge light inside 
~ 10" is dominant compared to the disk light. 

The construction proceeded as follows: The az- 
imuthally averaged luminosity profile from the ellipse fit 
to the observed I image (Fig [2]) was fitted by a two- 
component model including a steeper 7 - profile for the 
bulge and a natter 7 - profile for the disk. The contribu- 
tion of these profiles had to be equal at r « 10", and their 
sum had to fit the observed total profile. A good solution 
was found with 7 = 0.9 and a = 15" for the bulge com- 
ponent, 7 = 0.8 and a = 205" for the disk component, 
with a ratio of their total light B/D = 1.044. It is obvious 
that these two components have a rather similar struc- 
ture. Furthermore this decomposition is not unique and 
the astrophysical meaning of the labels 'bulge' and 'disk' 
for these components is not clear. 

Using these profiles we constructed the projected im- 
age of an elliptical prolate bulge with b/a ps 0.3 and its 
long axis along P. A. = 28° and an inclined disk (i = 35°, 
node line along P. A. = 113°) and coadded their projection 
as a twodimcnsional model image of M 94. The result was 
not a good description for M94, since the bulge compo- 
nent was too dominant to give the correct isophote twist 
at the correct radius. The same is true with a centrally 
flatter disk profile. Vice versa, if we weakened the bulge 
component or made its profile steeper, we could not sat- 
isfy above conditions b) and c). A number of different 
bulge axis ratios and bulge-disk combinations were tried. 
The most serious difficulty is that a bulge with a density 
profile as above has no sharp boundary, even though the 
bulge profile steepens to r~ 4 for r >> 15". Therefore no 
cusped isophotes are produced. A bulge does not have a 
natural outer radius like a bar at its corotation radius. Of 
course, it would be possible to use more complex compo- 
nents, e.g. a triaxial bulge with changing axis ratios or a 
bulge which is not elliptical or which has a sharp outer 
boundary. However, then too many arbitrary parameters 
would come into play and the result would not be very 
meaningful. 

We conclude that our first model is the most plausible: 
M 94 consists of a disk with a Hcrnquist like profile (7 w 
1) and a spheroidal bulge which can only be separated 
from the disk by its stellar velocity dispersion but not by 
its light profile or color. The deviations and twists of the 
isophotes are due to a superposed weak bar component. 
The bar ends inside the gas ring at its corotation radius. 

5. Dynamical model 

In this section we present a model for the kinematics in 
the inner region of M 94. We will try to explain the radial 
velocity curves of the stars and of the gas emission, espe- 
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daily of the 'overshooting' and 'counter-rotating' compo- 
nents (humps and dips, Fi g. p| ). 

As was shown in Sect. 3.7 the rotation curves of stars 



and gas are qualitatively similar. The velocity dispersion 
of the stars outside r = 20"is below 60 kms -1 , which is 
the observational limit. So we can assume that the kine- 
matics of the stars outside the bulge region is disk-like. 
For the bulge we have a projected stellar rotation veloc- 
ity v rot 50 kms _1 and a velocity dispersion a ~ 120 
kms -1 (Fig. Q). From the well known relations between e 
and v rot /(J we estimate a true flattening of e ss 0.4 for the 
bulge. 

5.1. Mass model 

We construct our mass model acc ording to the morpho- 



logical components from Sect. 4.1 



— The bulge and inner disk is approximated by a flat- 
tened 7 - model with 7 = 0.9, a = 15", and q = 0.6. 

— The bar component is approximated by a Ferrers bar: 
a = 20.8", 6 = 7.3", n = 2. 

— The contribution of the outer exponential disk for the 
mass model of the inner bar region is neglected. 

The mass-to-light-ratio £ was determined by fitting cir- 
cular orbits at r = 45" (inner ring) and was assumed to be 
equal for both components (no systematic color gradient 
was observed). 

The flattened 7 - profile for the inner disk and bulge 



Pg(p) 



with p 2 
10 10 M P 



M a 



2irqp~* (p + a) 4 -T' 



(3) 



z /q , scale-length a, and M — 2 x 



The mass of the bar can be calculated by integrating 
the bar surface brightness (Eq. [l]) over all space: 



M = 



27rCS a& / dr r(l - r 2 )™ 
Jo 

r(« + 1) 



T(n + 2)' 



(4) 
(5) 
(6) 



where the last integral was evaluated by using integral 
tables (e.g. Bronstein & Semendjajew, 1987, 1.1.3.4 #39). 

Using the values above one gets a bar mass of 9.8 x 
1O 8 M , which is about 14% of the mass of the flattened 
7 - model up to the end of the bar (r = 20" in the plane 
of the galaxy, 7.0 x 10 9 M Q ). From the surface photome- 
try calibration we obtain / = 7.8 inside that circle. With 
the color indices of de Vaucouleurs & Longo (1988) for 
the inner 35" aperture of M94, we obtain B = 10.6 inside 
of r = 20", and with a distance of 6.6 Mpc this leads to 



M/Lb = 1-8. This extraordinary low value was already 
obtained by Buta (1988); it is a consequence of the ex- 
tremely high surface brightness of M94. Fig. [i"2| shows the 
radial surface density of the different components of our 
model. 

The bar rotates rigidly with a constant pattern speed 
fib, which is determined in Sect. 



5.4 



Fig. 12. Surface brightness of 7 - model, Ferrers bar and their 
sum along the bar major axis in units of the characteristic 
length a. 



5.2. Gravitational potential 

To calculate the families of stable closed orbits one needs 
the potential of the galaxy model, from which one can de- 
termine the forces exerted on a test particle. The potential 
of the flattened 7 - model (Dehnen & Gerhard 1994) is 



* q {R,z) 

where 



1 



ip(rh)dT 



2 Jo (r + l)^^ 



(7) 



R 2 



t + 1 t + q 2 



and 

ip(rh) 



+ (3 - 7)m : 



(m + 1) 3 -t 
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This expression was numerically calculated using a 
routine by Walter Dehnen, whereas the potential of the 
Ferrers bar was calculated with the following expansion 
method described by Binney & Tremaine (1987): 



$(r, ip) 



oo k 

E E 



Vfc m Yfcm (arcsin(r), tp), 



(10) 



k=0 i 



-—k 



where the expansion constants Vkm are given by 

Vkra = -TT 2 Gg km / dip' 

Jo 

x / Z(r' : <p')Y k * m (zrcsm( r '),p')r'd r '. (11) 
Jo 

The functions Yk m are the well known spherical harmon- 
ics, G is Newton's constant of gravity and the gk m are 
constants. The variable r is normalized in such a way that 
the density distribution lies inside a circle of radius 1. 
R(<p') denotes the edge of the density distribution. The 
equations give the potential in a disk inside this circle. 
The expansion was truncated at k = 30 and m =8; such 
high values were necessary to approximate the potential 
of the elongated bar in M 94 accuratly enough (the axis 
ratio is 1:2.85). For rounder bars smaller values of k and 
I ml are sufficient. 



5.3. Resonances and families of closed orbits 

The families of closed orbits in barred potentials are deter- 
mined mainly by the loci of the Lindblad- and corotation- 
resonances (Contopoulos & Grosb0l, 1989). For weak bars 
linear theory can be applied. Then the resonances can be 
calculated in terms of two natural frequencies: 



1 rf$ 
~R~dR 



(12) 



is the circular frequency of a star on a circular orbit with 
radius Rq in the axisymmetric part of the potential ($0)- 
Small radial displacements from this circular orbit lead to 
radial oscillations with the epicycle frequency k , where 



dR 2 



ml 



(13) 



The frequencies fio and kq are calculated only from the 
axisymmetric part of the potential and are accurate only 
for weak bars. The exact positions of the resonances in 
strong bars have to be determined from the behaviour of 
the families of closed orbits (see next section). 

Azimuthal displacements do not change the circular 
orbit, so the azimuthal frequency is 0. The resonances arise 
between the forcing frequency seen by the star, m(f2o — 
Qb) (for barlike potentials m = 2), and the two natural 
frequencies kq and 0. 



If — Qb = 0, the Corotation resonance (CR) ap- 
pears, because the star moves as fast as the rotating bar. 
If m(r2 — ^b) = ±k , the star encounters successive crests 
of the potential at a frequency that coincides with that of 
the natural radial oscillations. An Inner Lindblad reso- 
nance (ILR) occurs if the star overtakes the rotating bar 
C+'-sign), an Outer Lindblad resonance (OLR) otherwise 
('-'-sign). 

A resonance- diagram can be plotted (Fig. |l3| ) and the 
loci of the resonances can be determined by drawing a 
horizontal line at The small oscillations in the curves 
near r = 40" are numerical artefacts: To calculate kq one 
has to evaluate the second derivative of <I>o- This axisym- 
metric part of the potential consists of the flattened 7 
- model and the axisymmetric part of the bar potential, 
which is given in a series. Even if the potential is approx- 
imated quite well, the second derivative still has oscilla- 
tions, which are visible in Fig. fL3[ 



Fig. 13. Resonance diagram for our model. The horizontal 
line Qt, = 2.5 shows the resonances, which are marked below. 
The bar ends at r = 20". The oscillations in k are due to the 
numerically calculated second derivative of the potential with 
a limited order of the expansion. 



The resonance diagram shows that because of the steep 
potential of M 94 only one inner Lindblad resonance can 
occur. In such a potential the closed orbit families in the 
corotating frame of reference have the following properties 
(see Contopoulos & Grosb0l, 1989 for a complete review): 
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The prograde xi-orbits are elongated along the bar and 
exist from the center outwards to nearly CR. These 
orbits should strongly sustain the stellar bar shape. 

The prograde X2-orbits are elongated perpendicular to 
the bar and exist in the energy range all the way from 
the center up to the energy of the ILR. 

The prograde a^-orbits have the same properties as the 
a^-orbits, but are unstable, so they can not be popu- 
lated with gas and are not considered in the following. 

The retrograde a^-orbits are not useful inside the CR, 
because they would lead to a retrograde streaming of 
gas, which is not observed. 

However, outside the CR we have only retrograde or- 
bits, which are nearly round. Inside the OLR they are 
slightly elongated perpendicular to the bar, outside the 
OLR along the bar. These orbits can be populated with 
gas, because a retrograde orbit outside CR has no ret- 
rograde streaming in the rest frame. 



u|j = v 2 — a 2 , v-rv z = 0, and = v 2 ot + er 2 , and the 
last term in the previous equation vanishes. Using 

<9$ v 2 - 

a* = if < 2 = °>- < 15 > 

one gets 

2 2 2 f d\np d\n<r 2 \ 

v ^ = v rot -° (^ + 7^. (16) 

where a is the isotropic stellar velocity dispersion, v rot the 
(observed) rotational stellar velocity, and v C i rc the circular 
velocity in the z = plane in the potential corresponding 
to the density distribution p(R,z). This equation is for- 
mally identically to that obtained from the Jeans equation 
in spherical symmetry, but is only valid in the z — plane. 

In this plane, using the density distribution of a Hcrn- 
quist sphere (Hernquist, 1990, corresponding to 7 = 1) 



5.4- Pattern speed 

The pattern speed f2b is the only free parameter in our 
model which can not be observed directly. Our aim was to 
determine fib using the structure of the orbits of our mass 
model. 

The stellar bar must be supported by xi-orbits. As 
these exist only inside CR, the pattern speed has to be so 
low that CR lies outside the bar (Sellwood & Wilkinson 
1992). Two other constraints on the pattern speed are the 
location of the gas ring at r ~ 45" and the gas-poor zone 
between 20" and 40". As was shown by Schwarz (1981) a 
rotating bar transports gas from the area between CR and 
OLR to or just beyond the OLR. To meet these two con- 
straints, the pattern speed in our model was determined 
at fif, = 2.5 in our units, corresponding to one revolution 
in 2.1 x 10 7 years or 290 kms _1 per kpc. In this way the 
OLR of the inner bar coincides with the ILR of the outer 
oval disk: the star formation ring at r « 45" (Gerin et al 
1991). 



5.5. Kinematics of the stars 

How do the observed stellar rotation velocity and the ve- 
locity dispersion in the bulge correspond to the observed 
mass distribution? The radial Jeans equation for an ax- 
isymmetric stationary system in the equatorial plane 
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(Binney, Tremame 1987, p. 197). Here R denotes the cylin- 
drical radius. Assuming an isotropic rotator, we have 



2nr (r + a) 3 
we obtain 

4- 



3a 
R + a 



(17) 



(18) 



neglecting the gradient of a 2 compared to the gradient of 
the density. The coefficient of a 2 in this equation goes to 
1 for R -> and to 5/2 for R -> a. 



The right hand side of Equ. (18) is obtained from ob- 
servational quantities of the stellar velocity field. We take 
the mean stellar rotation curve from Sect. 3.7 and depro- 
ject it using the optimal inclination of i = 30°. We assume 
a to be isotropic and spherically symmetric (no depro- 
jection) and calculate the mean velocity dispersion curve 
from the 4 independent stellar spectrograms (Fig. ||). The 
corresponding values oiv 2 irc in Eq. |l8|are shown in Fig. |IJ 
as circles. These are compared with the circular velocity 
curve from the axially symmetric part of the mass model 
(solid line in Fig. |l4j). Since the M/L ratio is a free param- 
eter here, the amplitude of the circular velocity curve was 
adjusted to the observed points. The important conclusion 
here is that observed points and the circular velocity have 
the same radial profile; the stellar kinematics confirm our 
mass model. 

For comparison we plot in Fig. fL4| b the inclination- 
corrected mean stellar rotation curve (+ symbols). This 
curve obviously has a shallower gradient and a lower peak 
rotational velocity, since it does not take into account the 
kinetic energy in a. We also plot the inclination corrected 
(i = 30°) gas rotation curve at P. A. — 118° (triangles), 
which is near to the line-of-nodes (113°) The maximum 
gas rotation velocity reflects the mass model fairly well. 
We see that the gas rotates faster than the stars, as ex- 
pected; however, the central gradient is again significantly 
shallower than that of the circular velocity curve. 
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clouds move on stable non-intersecting closed orbits (Sell- 
wood & Wilkinson, 1993). 

Figure |l5| shows the stable closed orbits in our model 
(inside CR only prograde orbits are drawn). The orbits 
are shown in the corotating frame of reference. All orbit 
integrations were done using a 4th order Runge-Kutta in- 
tegrator with variable time step and error prediction. This 
was mainly taken from Numerical Recipes in C (Press et 
al. 1992). 

In the very center of Figure |l5| one member of the X2- 
family is plotted, which shows the small extent of this 
orbit family. The ILR is at 5" which is small compared 
to the scale of the bar (20") and that of the figure. The 
xi-orbits are highly elongated along the bar and become 
nearly round near the CR. Between the CR and the OLR 
there exist no non- intersecting orbits, so gas is swept out 
of this zone and piles up at the OLR. Beyond the OLR 
there exist again non-intersecting closed orbits that can 
be populated with gas. These orbits are nearly round but 
slightly elongated along the bar. 



Fig. 14. a): The circular velocity, as calculated from the ob- 
served stellar rotation and velocity distribution by the Jeans 
equation (circles) , agrees well with the the circular velocity de- 

(solid and dotted curve, 



3.1 



rived from the mass model in Sect, 
resp.). b): Same as above, plus the stellar rotation curve (+ 
symbols) and the gas velocity (triangles) . 



5. 6. Kinematics of the gas 



The mass model of Sect. 5.1 was used to study the dy- 
namics of the optical line emitting gas. We neglect the 
self-gravity of this gas and consider it as test particles in 
the stellar potential. 

Gas has a high cross-section, so encounters between 
gas clouds normally have to be taken into account. These 
encounters lead to dissipation of energy and angular mo- 
mentum. As a first approximation we assume that the gas 



Fig. 15. Families of closed orbits in our model of M94. Res- 
onances are marked on the lower margin. While the center is 
populated by xi- and X2-orbits, in the region between CR and 
OLR no stable nonintersecting orbits can be found. Beyond 
OLR there exist nearly round retrograde orbits (in the coro- 
tating frame) which make up the ring. 



To obtain observable quantities from the orbits in Fig- 
ure [l5], one has to transform their velocities from the coro- 
tating frame to inertial coordinates and then project the 
orbits onto the sky. From the resulting velocity field we 
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calculate rotation curves at various position angles. Here 
we have to take into account that the velocity field is 
smoothed by atmospheric seeing (1.5" FWHM), by the fi- 
nite resolution of the spectrograph (120 km/sec FWHM), 
and by the finite slit width (4.3"). Therefore the gas ve- 
locity field has to be folded with the observed gas dis- 
tribution in the region of the galaxy contributing to the 
corresponding slit position. To obtain a model of the gas 
distribution, we consider an E-W cut of the H a distribu- 
tion (Fig. |]), this is shown in Fig. |l6|. We approximate the 
central emission peak by an axially symmetric distribution 
with a radial profile ~ (r + 2") , and the ring at radius 
r = 38"by a Gaussian profile (8" FWHM). 



Fig. 16. Emission gas density along a cut through the center 
of M 94. Three peaks in the gas density are visible. The inner 
peak corresponds to the center of the galaxy, while the outer 
ones at 45" are due to the star forming ring, which is shown in 
Fig. [|. The resonances in our model are marked on the lower 
margin. 



the circular velocity, which was used to calculate the mass 
of the model. However, between 5" and 10" the model 
curve shows a strong overshooting which arises from the 
contribution of the x\ orbits along the direction of their 
maximum streaming velocity. The observed rotation curve 
is less steep in the center and does not show such strong 
humps. This is similar for the adjacent slit orientations 
(c.f. Fig J). 

Fig. |l7b shows the model rotation curves along P. A. = 
28°, i.e. parallel to the bar. As could be expected the ro- 
tation curve along the bar is nearly flat since the bar axis 
is nearly perpendicular to the line of nodes. Our model 
does not show the observed humps at this slit position. 
However, the maximum and minimum velocities are rather 
high, since our broad slit always covers the innermost x\ 
orbits in their entirety. The slit smoothing gives the dis- 
played smooth curve, but this may depend on the precise 
distribution of gas that contributes to the observed veloc- 
ities. So the discrepancy in Fig. |i"7| b is not as serious as in 
Fig. |^a. 

Thus our kincmatical model is not able to reproduce 
the observed peculiarities of the gas rotation curves. From 
the model we expect humps in the curves especially along 
the minor axis of the bar (i.e. when the slit orientation 
is orthogonal to the bar) from the x\ orbits, and no 
humps on the major axis (spectrograms along the bar). 
Such a pattern is not observed, rather the distribution 
of humps and dips follows more circular patterns (spiral 
arms, Fig. O). 

As Fig. H shows, the gas rotation velocities are some- 
what higher than the stellar velocities, but resemble these 
much better than the curve in Fig. [Tjja. This is all the more 
surprising since in the inner 15" (the region of greatest 
discrepancy) the bulge is hot and the kinematics in agree- 
ment with the mass model (Fig. [Ti]) . This suggests that 
the H a gas in the inner 15" is both disturbed and still 
coupled to the stars. Cold molecular gas might well follow 
the model velocity field more closely, but corresponding 
high resolution CO data unfortunately are not available. 

Our kinematical model was constructed in such a way 
that the outer Lindblad resonance coincides with the gas 
ring at r w 45". Fig. [l8] shows the comparison of the ob- 
served gas velocity on the ring compared with our model 
calculations. 



Two example rotation curves were obtained by this 
procedure and are shown by the x symbols in Fig. 17 a,, b. 
The circles are the observed gas velocities along the cor- 
responding slit positions. The dotted curves give the min- 
imum and maximum velocities which occur in the consid- 
ered slit area before smoothing by seeing and finite reso- 
lution. 

Fig. [l7| a shows the model rotation curves along P. A. = 
118°, i.e. orthogonal to the bar. At radii greater ks 20" the 
agreement is quite good. This is because the orbits are 
nearly round there, so the velocities of the orbits reach 



6. Discussion 

6.1. NGC4736 - oval and barred 

The I image (Fig. [I]) and the ellipse fits (Fig. ^) show 
a clear nuclear bar in M94 (NGC4736), with axis ratio 
~ 0.3 and linear length ~ 30". The nuclear bar is im- 
plicit in the NIR data of Shaw et al (1993) and was also 
mentioned by Kormendy (1993). The galaxy NGC4736 is 
known as a prototypical oval disk galaxy. Such oval disks 
can be recognized both photometrically and kinematically 
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Fig. 17. Rotation curves along P.A. = 118° (top) and 
P. A. = 28° (bottom). The circles are the observed velocities, 
the crosses are the theoretical values. The dots give the min- 
imal and maximal velocity in every seeing element. For radii 
greater than m 20" the agreement is quite good. At the center 
the rise of the theoretical velocities is too steep to match the 
observations. 



(Kormendy 1982). Photometrically the brightness distri- 
bution consists of a series of distinct regions, which each 
have elliptical isophotes and a noticeable outer edge. The 
axial ratios and position angles of the isophotes in these 
nested oval regions are different, implying that at most one 
of them can be axisymmetric if they are all coplanar. Kor- 
mendy (1982) argues that warps are rare in edge-on spirals 
at the high surface brightnesses at which the oval struc- 



Fig. 18. Measured velocities of gas in the ring of M94 with 
respect to the center, as a function of posit ion angle. The 
asterisks represent our measurements in Sect. |3.8| . The other 
symbols indicate different emission lines from van der Kruit 
(1974). The dotted lines give the theoretical predictions based 
on closed orbits in the ring with approximate radii of 42", 47" 
and 51". 

tures occur (< 25-B-mag arcsec~ 2 ) and that therefore the 
observed features must be truly non-axisymmetric. Kine- 
matically, oval disks signify their presence by several fea- 
tures in their velocity fields which are also known from 
barred or triaxial galaxies: The kinematic major axis is 
different from the photometric major axis, it is not per- 
pendicular to the kinematic minor axis, and its position 
angle may change with radius. These signatures arise from 
the projection of quasi-elliptical stellar orbits or gaseous 
streamlines orientated at oblique angles with respect to 
the line-of-sight and line-of-nodes. 

The B-band surface brightness profile of M 94 was mea- 
sured by Boroson (1981); it appears to show several char- 
acteristic radii: 40" and 120". Our /, J, K brightness pro- 
files confirm the shelves at 40" (1.3 kpc) and 120", al- 
though the feature at 40" is less pronounced. The reason 
for this is probably the blue color of the star forming ring 
at r = 45". Numerous large HII regions are located on 
this inner ring which modify the blue surface brightness 
profile (Fig. ||). The inner ring is also seen in HI (at 40- 
50", Bosma et al 1977, Mulder & van Driel 1994), in CO 
(Garman & Young 1986, Gerin, Casoli & Combes 1991), 
and in H a (Fig. f§ van der Kruit 1974, Buta 1988). The 
outer disk of the galaxy is not obviously barred, but the 
run with radius of ellipticity and position angle suggests 
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an oval shape (Fig. |^). M94 is well known for its faint 
outer ring at ~ 330" seen in HI and the stellar distribu- 
tion (e.g. Plate I in Lindblad 1961 or Fig. lb in Mulder 
& van Driel 1993). From our kinematic and photometric 
study we think it is most likely that the inclination of M 94 
is i ~ 35°, so that the most axisymmetric region is that 
portion of the disk between 50" (1.7 kpc) and 100" (3.3 
kpc, Sect. ^2). 

The kinematics of the HI and H a gas at medium and 
large radii show significant departures from axial symme- 
try. A recent model of the outer oval disk is presented by 
Gerin et al (1991), who assume an oval shape of the po- 
tential between 5.5 kpc and 13 kpc and associate the two 
rings with the inner and outer Lindblad resonances (ILR 
and OLR) of this rotating oval. In Sect. 3.3 we concluded, 
that the disk of M 94 is oval at least between 120" (4 kpc) 
and 220" (7.5 kpc); the inner radius is compatible to the 
values of Gerin et al, our outer value is lower. However, 
since the ellipticity of the disk still increases between 220 
and 400", the oval disk might well be larger. 

The modelling presented in this paper, based on opti- 
cal CCD and NIR observations of the central bulge and 
bar, suggests that the most likely value for the pattern 
speed is such that the OLR of the bar coincides with the 
position of the inner HI and H a ring at 45". Combining 
this result with the simulations of Gerin et al (1991) it 
thus appears that M 94 is a galaxy with an inner bar em- 
bedded in an outer oval disk, such that the OLR of the 
inner component is located approximately at the ILR of 
the outer component. Such bar- within-oval galaxies share 
many of the characteristics of the bar-within-bar galaxies 
discussed, e.g., by Pfenniger & Norman (1990) and Friedli 
& Martinet (1993). In the following, we discuss some as- 
pects of this picture for M 94 in more detail. 

6.2. Nuclear bulge and nuclear bar? 

Does NGC 4736 have a nuclear bulge? On the basis of ear- 
lier photometry and kinematic measurements, Kormendy 
(1993) concluded that M94 is one of a class of galaxies 
that appear to contain central bulges which are really only 
disks. The evidence he presented in favour of this view was 
the following: (i) The central component with the r 1 / 4 
brightness profile ('bulge') coincides spatially with spiral 
structure and a nuclear bar. These features are charac- 
teristic for disks; thus the r 1 / 4 profile cannot per se be 
used as evidence for a truely three-dimensional bulge, (ii) 
Based on data by Pellet & Simien (1982), the apparent 
v/a ~ 0.77 of NGC 4736 is well above that expected 
for an oblate isotropic rotator of its flattening (e ~ 0.1, 
(v/a) h\ — 0.35). Thus the central kinematics appear disk- 
like rather than bulge-like. 

With the new data presented in this paper we can 
now investigate this question in more detail. The stellar 
velocity dispersion profile (Fig. [sj) shows a clear transi- 
tion from an outer cold disk to an inner hot component. 



Because we are observing M94 nearly from above, the 
measured dispersion mostly characterizes the vertical mo- 
tions in this galaxy. Outside 15" the observed dispersion is 
consistent with the instrumental resolution (~ 60kms~ 1 ) 
and at those radii the light may be dominated by a cold 
disk component. In this region the circular velocity is 
v c ~ 200kms _1 . Inside 15", the projected velocity dis- 
persion rises to 120 kms -1 , while simultaneously the mea- 
sured rotation velocity linearly decreases inwards to zero. 

From the (admittedly noisy) line-of-sight-velocity dis- 
tribution ('LOSVD') shown in Fig. |, we see that in fact 
the rotation velocity at 10" is still dominated by the 
narrow peak corresponding to the rapidly rotating disk. 
However, the bulge has a much broader, more symmet- 
ric LOSVD with a mean rotation velocity of ~ 50 kms" 1 . 
This value reduces the intrinsic bulge V max /a to a value 
of ~ 0.8 (deprojected), which is more in accord with its 
intrinsic flattening of q w 0.6 (deprojected with e = 0.1 
and i = 35°). 

If indeed the entire central light came from a constant- 
thickness disk, then the observed surface brightness pro- 
file, which is consistent with a projected Hernquist model, 
would imply a surface brightness oc R , and hence a ver- 
tical velocity dispersion cx Rr 1 ! 2 . Given our limited res- 
olution, the predicted rise in the vertical dispersion by a 
factor ~ 2 between 15" and 3" is consistent with the mea- 
sured value of the central dispersion. However, the shape 
of the velocity dispersion profile, in particular the clear de- 
cay at around 15", argues in favour of a three-dimensional 
central bulge component, supporting the analysis of the 
LOSVD above. 

We next discuss the nature of the central bar. From the 
photometry we infer that the bar comprises only ~ 15% of 
the light within 20". The light distribution is significantly 
more axisymmetric in the region r < 10" and r > 20" 
than near r = 15" where the bar is strongest in the sur- 
face brightness profile (Fig. |^). The new surface photom- 
etry shows that the bar signature disappears very quickly 
for r > 15". The brightness distribution cannot be fitted 
with a single ellipsoidal component as the isophotes are far 
too pointed near the end of the bar (04 profile in Fig. ||) . 
A strong gradient in the ellipticity and boxyness parame- 
ter 04 would be required in a one-component model. This 
rules out a single triaxial bulge with nearly constant axis 
ratio as an explanation of the data. However, it also ar- 
gues against a single disk-like component, for this would 
have to be almost completely axisymmetric in the inner 
10" and develop a separate thin (b/a ~ 0.3) bar compo- 
nent near r = 15", a morphology which is not usually seen 
in bulge-less, barred galaxies. 

Because the projected length of the bar (15"), the size 
of the bulge, and the scale of the linearly rising part of the 
rotation curve (Fig. ^) all coincide, we can not use the ve- 
locity dispersion data to determine the vertical extent of 
the bar. Indeed, the fact that these scales approximately 
agree might suggest that the bar and the bulge could after 
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all be one and the same component. However, the bright- 
ness profile of the bar is much shallower than that of the 
bulge, so that in this interpretation the quadrupole com- 
ponent would have to be strongly concentrated towards 
the outer parts of the bulge. Moreover, along its minor 
axis direction, the bar is significantly smaller (» 7") than 
the bulge (15"). The former property appears consistent 
with the properties of peanut/box-shaped bulges, which - 
compared to an ellipsoidal model - have more excess flux 
near their edges than in their central parts (Shaw 1993). 
The latter property can not be compared with these ob- 
jects because peanut bulges are usually observed edge-on. 

Taken together, the photometry and stellar kinematics 
of M 94 therefore appear to show that this galaxy contains 
a small central bulge, a centrally concentrated disk, and a 
weak 30" diameter bar with shallower profile than either 
of the former two components. The bulge may be super- 
posed on the central disk and bar, or, perhaps more likely, 
because of the similar surface brightness profiles of bulge 
and disk, simply be the inner, more spherical part of the 
disk. 

6. 3. Pattern speed of the central bar 

The model discussed in the previous Section assumed that 
the bar rotates at a pattern speed f2 p = 290kms~ 1 /kpc, 
such that corotation lies at 25". This places the inner ring 
of HI and H a gas at approximately the OLR of the rotat- 
ing potential; as shown by Schwarz (1981) gas between CR 
and OLR is given angular momentum by the bar, moves 
outwards and collects in a ring of material near the OLR. 
This zone interior to the OLR is characterized by a weak 
spiral pattern, which can best be seen in the morphology 
of the dust lanes (Fig. ||). The corotation radius then cor- 
responds approximately to the end of the bar in the CCD 
image (Fig.g and Fig. |o|), and to the region of lowered 
cold gas density visible in the CO observations of Gerin 
et al (1991). The depression near corotation is extremely 
pronounced in the H a emission distribution (Fig. || and 
Fig. 0). 

Although a model with this value of f2 p thus explains 
a number of features observed in M94, it failed in the 
sense that the measured H a gas kinematics in the central 
15" are not fitted by those predicted if the gas moved on 
closed orbit streamlines in the barred potential. However, 
this may be more a problem to do with the H a gas: in- 
deed for no reasonable value of the pattern speed was it 
possible to obtain a satisfactory fit to the H a kinematics. 
In any case it is unlikely that the bar could rotate faster, 
for then corotation would be placed inside the bar. There 
would be no orbits to support the outer bar, yet the stellar 
population of the bar appears old (> 3 Gyr). 

Could the bar rotate more slowly? This at least would 
go in the sense of somewhat easing the discrepancy of the 
predicted and the H a kinematics. In this case the inner 
ring at 45" might be at ILR rather than at OLR, and 



then the bar's pattern speed would be equal to the pat- 
tern speed of the outer oval in the Gerin et al model. If 
this were the case, the bar would end inside its ILR. How- 
ever, such a model could only be in dynamical equilibrium 
if the principal axes of bar and oval were either aligned or 
perpendicular to each other (Louis & Gerhard 1988), nei- 
ther of which possibility is in accord with the observed 
misalignment angle between bar and adjacent disk: From 
P.A. bar = 28° and P.A. ovai = 90° we get A « 62° in the 
sky or ~ 58.5° in the plane of the galaxy for i = 35°. If, 
on the other hand, the two components are oriented at 
an oblique angle, then torques between them will have a 
strong effect particularly for low pattern speeds. Thus we 
conclude that the value of tt p = 290kms~ 1 /kpc that we 
have used is the most likely one. 

Models with two non-axisymmetric components rotat- 
ing at different pattern speeds were discussed and gener- 
ated by Pfenniger & Norman (1990) and Friedli & Mar- 
tinet (1993). These authors concluded that in order to 
avoid too large resonant and chaotic zones in the orbit dis- 
tribution it was favourable to place the ILR of the outer, 
more slowly rotating component onto the CR resonance 
of the inner, more rapidly rotating bar. Indeed, in the nu- 
merical (N-body and SPH) simulations of Friedli & Mar- 
tinet (1993) two-bar systems with the corresponding ratio 
of pattern speeds formed. In the formation of the inner 
secondary bar, dissipation is found to be critical in assem- 
bling sufficient amounts of gas onto the X2-orbits in the 
primary bar potential. 

However, the two-bar phase in these models only lasts 
for about five rotations of the faster, inner bar, after which 
this dissolves. Therefore, if this is to be a viable sce- 
nario for M 94, the inner bar must have been triggered 
very recently, e.g., by inf ailing molecular gas, and it must 
currently end somewhat outside its ILR. This appears 
difficult in the dominant field of a nearly axisymmetric 
component. Moreover, the stars themselves in the M94 
bar are several times older than the estimated dissolution 
timescale. Thus one would require that during their bar 
triggering no young stars, formed in the molecular clouds, 
should have been mixed into the stellar bar, which ap- 
pears unlikely. Alternatively, if the inner bar indeed ro- 
tated faster, such that the ring is at OLR, then the inter- 
action with the outer oval would be weaker because of the 
adiabatic invariants, and the lifetime of this configuration 
could be longer. 

6.4- Dynamics of the inner H a - emitting gas 

Warm optical emission gas can be detected in M 94 mainly 
in the inner star forming ring (r f» 45") and in the central 
bar region. The dynamics of the ring was studied in detail 
by van der Kruit (1974, 1976) and Buta (1988). The veloc- 
ity field of the innermost region is fairly complex. There 
are several aspects: 
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(1) The H a streaming velocities measured in the cen- 
tral 15" of M94 are very similar to the stellar streaming 
velocities in the bulge (Fig. ||). These correspond to a hot 
stellar system, as the deprojected rotation is v 2 ~ a 2 /2 



at 10" (Sect. 5.5). Thus the H a streaming velocities are 



significantly less than the circular velocity estimated from 
the Jeans equation. However, the gas rotates faster than 
the stars for r > 20", fitting to the HI rotation curves of 
Mulder & van Driel (1993). 

(2) The gas velocity field does not reflect the high 
streaming velocities expected from closed orbit stream- 
lines in the barred potential. In particular, the expected 
steep rise of the measured velocities in the central 5" for 
slit positions orthogonal to the bar is not reproduced. 

(3) The warm gas velocity field shows, superposed on 
the general near-circular rotation field, a number of humps 
and dips which seem to correlate with the weak dust spiral 
pattern observed in the region between central bar and 
ring (Fig. |). 

We now discuss several possible explanations for the 
observed low gas velocities in the very center. 

(a) The gas velocities are intrinsically higher, but we 
do not observe the high values because of low spatial res- 
olution (seeing and broad slit) combined with the finite 
spectral resolution. These effects were taken into account 
for the construction of the rotation curves in Fig. |l7], so 
this possibility can be rejected. 

(b) The gas in the region of the central bar moves 
intrinsically similar as the stars. This might be plausible 
if the observed gas originated from mass loss from the old 
stars, and still shared their velocity field. 

(c) Walker et al (1988) published 2.0 — 2.5/im spectra 
from the central 8" region of M 94. They showed that the 
B 1 emission is very weak there compared to typical star- 
burst galaxies like M82 or NGC253. Therefore no strong 
star formation is presently going on in the center, con- 
sistent with the not very strong H a emission from the 
central region. However, from the still relative strong CO 
absorption bands of young supergiants Walker et al (1988) 
concluded that a strong starburst did occur in the central 
region of M 94 some 10 7 years ago. A burst of star forma- 
tion produces wind bubbles and supernova shocks unre- 
lated to the potential of the galaxy. Presently, the nuclear 
region of M 94 contains a mixture of warm gas, warm dust, 
and molecular clouds (Smith et al 1991). The warm gas 
might not yet have settled to cloud orbits in the plane, 
making its present partly chaotic velocity field plausible. 

It will be interesting to compare future high-resolution 
molecular gas kinematics with the model streamlines in 
Sect. 0. 



7. Conclusions 

(1) We have reported CCD imaging in / and H a , focal re- 
ducer imaging in V, I, and NIR J, K imaging of the bright 
Sab galaxy M94 (NGC4736). This galaxy is known as 



an oval disk galaxy. Two rings, an outer ring observed at 
330"and an inner ring at 45", are thought to be located 
near the outer and inner Lindblad resonances (OLR and 
ILR) of this rotating oval (Gerin et al 1991). Our study 
of the surface brightness profiles and isophote shapes con- 
firms the oval disk structure of M94 between r = 120" 
and r = 220". Basic parameters of M 94 and the most im- 
portant numerical results from our work are collected in 
Table 1. 

(2) The central surface brightness distribution of M 94 
is characterized by a power-law profile. It is not possible to 
discriminate bulge and disk components from the photom- 
etry alone. We have constructed a morphological model of 
the inner regions in terms of an axisymmetric model with 
a 7 = 0.9 cusp and scale-length a — 15" (Dehnen & Ger- 
hard 1994). Outside r — 60" a nearly exponential disk 
with scale-length f3 w 60" must be added to this model. 

(3) Inside the inner ring the position angles of the 
isophotes rotate by nearly 90°. The isophotes show cusped 
deviations from ellipses in a range of radii just interior 
to r = 15". After subtracting the axisymmetric model, 
the residuals show a nuclear bar, with projected radius 
r = 15", axial ratio b/a = 0.3 and P. A. = 28°. The outer 
oval disk has P. A. — 90°; the deprojected misalignment 
angle between these two components in the plane of the 
galaxy is A = 58.5°. The luminosity of the inner bar is 
5.9-lO 8 Z/0 (in B), comprising about 14% of the total light 
inside 20". 

(4) Thus M94 has a bar-within-bar morphology. We 
suggest that the corotation radius of the inner bar is just 
outside its deprojected major axis length (20"= 0.64 kpc). 
This places the inner gas ring (at 45"), which is believed 
to fall near the ILR of the outer oval, at the OLR of the 
central bar. The zone near corotation (CR at 26") is char- 
acterized by a marked depression in the H a and HI gas 
densities, consistent with the absence of non-intersecting 
closed orbits in this region. Between CR and OLR of the 
central bar, spiral arms are detectable, especially well by 
their thin dust lanes in, e.g., I — K color index images. 

(5) We have presented absorption line spectroscopy in 
4 longslit orientations in the inner 60" of M 94. The stellar 
velocity dispersion profile reveals the existence of a cen- 
tral bulge of 15" radius, with a — 120kms _1 . Outside the 
bulge we observe a — 60kms~ 1 in the disk region, which 
is the spectral resolution. The stellar rotation curves are 
smooth and are consistent with an axially symmetric ve- 
locity field with inclination i = 30° and a line of nodes at 
P. A. = 113°. This suggests that the effects of the nuclear 
bar on the stellar velocity field are weak, in the bulge re- 
gion because this is a hot stellar system, and outside the 
bar because there the potential becomes rapidly axisym- 
metric. The combination of stellar rotation and velocity 
dispersion is consistent with the photometric model in the 
bulge and adjacent disk regions in the sense that the Jeans 
equation is satisfied for M/Lb ~ 1.8. 
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(6) Futher evidence that the bulge and disk compo- 
nents in M 94 are distinct comes from the line-of-sight 
velocity profiles at positions along the minor axis but 
beyond the extent of the nuclear bar, which suggest a 
rapidly rotating, cold component superposed onto a more 
slowly rotating, hot component. From these profiles we es- 
timate that the deprojected mean rotation of the bulge is 
v w lOOkms -1 , and that v/a = 0.8, roughly consistent 
with an intrinsic bulge axis ratio c/a = 0.6. 

(7) By contrast, the kinematic data do not allow us 
to determine whether the bulge and bar components are 
distinct. The result that the radius of the bulge and the 
major axis length of the bar are similar might suggest that, 
in fact, bulge and bar are one and the same component, 
perhaps similar to that inferred in the inner Galaxy. How- 
ever, the bar is significantly smaller along its minor axis 
direction than the bulge, so that a simple near-ellipsoidal 
component is not a good model for the data. 

(8) We have obtained emission line spectroscopy in 8 
longslit orientations of the inner 100" of M94, and thus 
obtained new results on the gas kinematics in the region 
of the inner bar. The streaming velocity field of the 10 4 K 
gas is similar to that of the stars for r < 20". The gas 
rotates faster than the stars and fits the HI rotation field 
(Mulder and van Driel 1993) at larger radii, consistent 
with circular orbits. 

(9) Orbit calculations in the gravitational potential 
derived from the photometry of M94 predict large non- 
circular motions for cold gas in equilibrium in the bulge 
region. The resulting kinematics do not reproduce the ob- 
served central H a kinematics. For the future it would be 
extremely interesting to compare high-resolution CO gas 
velocity observations with the model predictions. 

(10) The result that the H a streaming velocities in the 
inner 15" are less than required for gravitational equilib- 
rium suggests that the warm ionized gas has significant 
random motions. This could be because this gas is still 
associated with the stars through recent mass loss, or, 
perhaps more likely, that it is beeing heated by shocks 
and localized gas flows. Indeed, while Walker et al (1988) 
found no B 1 in the central 8" in M94, i.e. no signs for 
an on-going starburst, they concluded from the strengths 
of the CO bands that a strong starburst must have oc- 
curred several 10 7 yr ago. Direct evidence for some still 
ongoing star formation in the inner regions of M 94 also 
comes from the observed H a emission. 
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Table 1. Parameters of M94 



Item 



Quantity 



Reference 



General 
Parameters 



Morphology 
on Sky 



Inner Disk+Bulge 
Surface Photometry 



Bulge Model 



Central Bar Model 



Stellar Kinematics 



type 

distance 

scale 

Vsys 

mass 
B t 

bulge 

central bar 
inner ring 
oval disk 
outer ring 

brightness profile 

ellipticity 

position angle 

inclination 

mass of inner disk 

flattening 
mass 
v/a 
M/L B 

axis lengths (deprojected) 
mass 

pattern speed 
ILR inner Lindblad res. 
CR corotation radius 
OLR outer Lindblad res. 

v rot (deprojected) 
velocity dispersion 
inclination 
line of nodes 



(R)Sab(s) resp. (R)SA(r)ab II 

6.6 Mpc 

I" = 32pc 

310 km/s 

5.1 x 1O 1O M 

8.96 mag 

r = 15 - 20" 

a = 15", b = 7", P.A. = 28° 
r = 45", P.A. = 127°, e = 0.27 
r = 120 - 220", P.A. « 95°, s = 0.23 
r « 330", P.A. « 150° 

7 = 0.9, a = 15" 
e « 0.18 
P.A. « 110° 

i w 35° 



7.0 x 10 9 M m 



g = 0.4 
2.0 x 10 9 M, 
0.8 
1.8 







a = 21", b = 7" 
9.8 x 1O 8 M 
290 km/s/kpc 
5" = 160pc 
26" = 830pc 
40" = 1280pc 

230km/sec max. at 15" 
120 km/sec 
i = 30° 
P.A. = 113° 



RSA resp. RC2 



Buta 1988 

Mulder & van Driel 1993 
Buta 1988 

Sect. 3.7 
Sect. 3.1 

van der Kruit 1976 
Sect. 3.2 
Buta 1988 



Sect. 
Sect. 
Sect. 
Sect. 
Sect. 

Sect. 
Sect. 
Sect. 
Sect. 



3.2 
3.4 
3.4 
3.4 
5.1 

3.7 
5.1 
3.7 
5.1 



Sect. 5.1 

Sect. 5.1 

Sect. 5.4 

Sect. 5.4 

Sect. 5.4 

Sect. 5.4 

Sect. 5.2 
Sect. 3.7 
Sect. 3.7 
Sect. 3.7 



Gas Kinematics 

HI 

Molecular Gas 



Vrot (deprojected) 
mass 

Vrot 

inclination 
line of nodes 

Hi mass 

Vrot 



250km/sec max. at 15-20" 
1O 9 M 

130 - 170fcm/s 
i w 40° 
P.A. « 115° 

5.1 x 10 s M Q 
200 km/s at 60" 



Sect. 3.8 

Mulder & van Driel, 1993 
Mulder & van Driel, 1993 
Mulder & van Driel, 1993 
Mulder & van Driel, 1993 

Gerin et al 1991 
Gerin et al 1991 



